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COMPLETENESS AND CONVERGENCE 


C. D. U. 517.96 


IN THE FINITE 


ELEMENT METHOD (*) 


RESUMO 


Trata-se de um estudo sobre a convergência para a 
solução exacta de sucessões de soluções aproximadas 
geradas por sistemas de elementos finitos com dimen- 
sões decrescentes. Justifica-se um critério geral de com- 
pletidade, e demonstra-se que um tal critério o é também 
de convergência se as derivadas de segunda ordem dos 
deslocamentos correspondentes às mesmas soluções apro- 
ximadas permanecerem limitadas dentro de cada ele- 
mento, à medida que as dimensões respectivas tenderem 
para zero. As conclusões são válidas para todas as teo- 
rias estruturais lineares continuas, isto é, para a Elas- 
ticidade bi e tridimensional, assim como para vigas, 
cascas e placas. No caso das teorias simplificadas, isto é, 
das teorias em que se despreza a deformação por esforço 
transverso, as terceiras derivadas do deslocamento trans- 
versal devem também permanecer limitadas dentro de 
cada elemento, para que a convergência seja conseguida. 


SECTION I 
INTRODUCTION 


The problem of convergence to the exact 
solution of a sequence of approximate solutions 
generated by patterns of finite elements decrea- 
sing in size is the main concern of this paper. 


por EDUARDO R. DE ARANTES E OLIVEIRA 


Professor of the Instituto Superior Técnico, Head of the 
Applied Mathematics Division, Buildings and Bridges 
Department of the Labcratório Nacional de Engenharia 
Civil, Lisboa. 


SYNOPSIS 


À study is made of the convergence to the exact sotu- 
tion of seguences of apbroximate solutions generated by 
elements with decreasing size. A general completeness 
criterion is justified which is shown to become also a 
convergence criterion if the second order derivative of 
the displacements corresponding to the sucessive ap- 
proximate solutions remain bounded within each element, 
as the size decreases indefinitely. The conclusions are 
valid for all the different linear continuous structural 
theories, i.e, for two and three-dimensional Elasticity, 
as well as for beams, shells and plates, If the simplificd 
theories which result from neglecting the transverse 
shear deformation are considered, the third derivatives 
of the transverse displacement must also remain bounded 
within each element, in order that convergence can be 
obtained. 


In the Ritz method (Reference 1), completeness 
is a sufficient condition for convergence. It is 
important to notice that the finite element method 
is a particularization of Ritz method only if 
compatibility between elements is achieved 
(Reference 2). Convergence is however still 
possible if compatibility is violated. Completeness 


(*) First published in the Proceedings of the 2nd Conference on Matrix Methods in Structural Analysis, Dayton, 


Ohio, 1968. 


TECNICA Nº 403 


109 


and convergence criteria are justified in the 
present paper which are valid even if compati- 
lity between elements is not obtained. The 
different kinds of structures are simultaneously 
considered but a special Section will be dedicated 
to the simplified theories which result from 
neglecting the transverse shear deformation. 

A general discussion of the subject was already 
presented in a former paper by the author 
(Reference 2). The finite element method, how- 
ever, was considered as a general mathematical 
technique, while the present paper is only 
concerned with elastic problems. Such particular- 
ization makes the discussion less abstract and 
thus easier to follow. 

Some other papers (References 3, 4, and 5) 
have been written which examine convergence 
in connection with particular kinds of elements. 
No general convergence criteria have however 
been presented or justified in those papers. 


SECTION II 
A SYNTHETIC FORMULA TION 
OF ELASTICITY 


Elastic theories involve three kinds of magni- 
tudes, stresses, strains and displacements, whose 
vectors will be denoted by s, € and u. 

Such magnitudes are defined on a domain, D, 
corresponding to the body, and related by three 
kinds of field equations which can be symbolized 
as follows: 


a) Equilibrium equations: 

Es=X (1) 
b) Strain-displacement relations : 

e=Du (2) 
c) Stress-strain relations: 

s=Hoe (3) 


E and D are first order differential operators, 
X is the vector of the body force density compo- 
nents, H is a symmetric positive definite matrix. 

Equations 1, 2 and 3 are valid on D. On the 
boundary, B, the equilibrium equations become 


Ns=p (4) 


N is a matrix whose elements depend on the 
orientation of the normal vector at each bound- 
ary point. p is the vector of the tractions 
applied to the boundary. 

The elastic analysis reduces to finding the 
solution of the system of field Equations 1, 2 and 
3 which satisfies certain boundary conditions. 
The simplest and most important types of 
boundary conditions can be expressed directly 
in terms of displacements or tractions applied 
to the boundary. Let B, and B, denote the 
portions of the boundary where tractions or 
displacements are respectively prescribed. The 
boundary conditions are analytically expressed by 


(5) 
u=u in B, (6) 

The stresses and strains can be eliminated an 
the problem can be formulated in terms of 


displacements only. Indeed, introducing Equations 
2 and 3 into 1, there results 


Ku=X (7) 
in which K is the second-order differential 
operator 

K=EHD (8) 
The boundary conditions are also easily 


expressible in terms of displacements or displace- 
ment derivatives. 

Any set of fields, 8, e, u, defines an elastic 
field. 

An elastic field is said to be compatible if 
Equations 2 and 3 and boundary conditions 
Equation 6 are obeyed. An elastic field is said to 
be equilibrated, with respect to a certain set 
of external forces, if it respects Equations, 1, 3 and 
4 and boundary conditions Equation 5. The exact 
solution is the elastic field which is simulta- 
neously a compatible and an equilibrated one. 

Assume now that u and $s are such that the 
derivatives involved in D and E exist everywhere 
in D. Then, the following relation holds 


e T T 
s (Du)dD= | (Es) udD+ | (Ns) udB 
Jem nam + [mm neo 


D D 
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In this relation, vectors S and U are not neces- 
sarily related by the stress-strain relations, 
Equation 3. 

The theorem of the minimum total potential 
energy, which is a consequence of Equation 9, 
states that the exact solution minimizes the 
functional 


TT = 
F=U- [X vdD— fp dB (o) 
D B 


i.e. the total potential energy, in any class of 
compatible elastic fields which contains the exact 
solution. 

The first term in EF, 


A eis: 
U=[waD=— [e HedD= 


D D 


1 T 
= [ (Du) H(Du)dD (11) 
D 


is the strain energy. 


The upper dash on X and p indicates prescri- 
bed magnitudes. 

The formulation which has been presented is 
quite general as it is valide not only for linear 
two-and three-dimensional Elasticity, but also for 
linear theories of plates, shells and beams. 

In the case of a plate, for instance, vector U 
contains the transverse displacement and two 
rotations, vector € contains the curvatures and 
the transverse shear strains, vector S contains the 
bending ande twisting moments and the trans- 
verse shearing forces. 

A very frequent simplification in the analysis 
of plates, shells and beams consists in neglec- 
ting the transverse shear deformation. This 
makes it possible to reduce the number of un- 
knowns to one (the normal displacement) in the 
theory of plates, and three (the normal displace- 
ment and the tangential displacements) in the 
theory of thin shells. 
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Such simplified theories present the same 
scheme as the initial ones. Operators E and D 
are now however of the second order and oper- 
ator K of the fourth order. The compatibility 
conditions involve the continuity of both the 
displacements and their first derivatives. As the 
rotations are expressed in terms of such magni- 
tudes, the compatibility conditions are essentially 
the same, both in the simplified theories and in 
the corresponding initial ones where the trans- 
verse shear deformation is not neglected. 

In the following discussion, the initial theories 
will first be considered. Operators E and D will 
thus be assumed of the first order, and continuity 
of displacements will be the only condition for 
compatibility Then, at the end of the paper, the 
results established for the initial theories will be 
used in the discussion of the simplied theories. 


SECTION IN 
THE FINIT ELEMENT METHOD 


The finite element method is a general discre- 
tization technique. In this method, domain D is 
assumed to be decomposed into a finite number 
of subdomains D“ and families of fields are con- 
sidered which have different analytical expres- 
sions within each subdomain. 

A finite element is a closed subdomain, D”, 
together with the family of elastic fields which 
are allowed to occur within it. Such family is a 
linear combination with coefficients q, of a finite 
number of unit modes, so that each field of the 
family corresponds to ascribing particular values 
to the parameters di: 

The values of the displacements at a certain 
number of points on the boundary of the element, 
called nodes or nodal points, are, as a rule, 
chosen as parameters. The allowed fields need 
not be introduced however by giving the expres- 
sion of the displacements directly in terms of 
their own nodal values. They can indeed be given 
in terms of equal number of arbitrary parameters 
which in turn can be expressed in terms of those 
nodal valures. 

The type of an element refers to its general 
shape, nodal point specification and to the 
allowed fields, analytically defined by expressing 
the displacement vector uí in terms of the para- 
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meters (generalized displacements) and the coor- 
dinates : 


e 


v=f(x,x ++ )Q (12) 


qó is the vector of the nodal displacements gs. 


The elements É; of fº are supposed to be 
continuous in the closed subdomain occupied by 
the finite element. The unit modes are defined 
by the columns of f”. 

Each displacement component is assumed to 
depend only on its own nodal values. Thus, if 
q; corresponds to the component u; , allthe mag- 
nitudes ta; for which k =, will vanish. 

In order that Equation 12 can be homogeneous, 


q; must be amenable to a nondimensional form 


mdf 58 
Jedi (o + +++) (13) 


in which [É is a typical dimension of the element, 
for instance its maximum diameter. vi is a 
function which does not depend on the absolute 
dimensions of the element. 

The different finite elements are made com- 
patible through the specification of reduced com- 
patibility conditions. These require that the values 
of the displacements be the same at coincident 
nodes of adjacent elements and equal the pre- 
scribed ones at the nodes located on B,. For the 
sake of simplicity, such prescribed displacements 
are assumed to vanish. 

A point of the domain is said to be a node 
of the system if it is a node for on or more 
elements. 

Let q, be the vector of the values of the 
displacements at every node of the system but 
those which are located on B,. The reduced 
compatibility conditions can be expressed by 
writing for each element 

q=Ta4, (14) 
where matrix Tº depends on the topology of 
the system. 
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Equations 12 and 14 show that the knowledge 
of q, is enough for the definition of the field 
within every element of the system. 

Introducing 14 in 12, there results 


uº = Fº q (15) 


in which 


E = 107" (16) 


Equation 15 is valid in D*. The whole set of 
equations (15) (one for each element) piecewise 
defines a family C, of elastic fields each of 
which corresponds to a certain vector q, . Index 
n refers to a certain degree of subdivision of D 
into subdomains. 

The reduced compatibility conditions are gen- 
erally not sufficient to make the displacement 
components continuous across the element bound- 
aries. This depends on the type of the element. 

If the type is such that the reduced compatibility 
conditions are sufficient to ensure continuity of 
the displacements across the element boundaries, 
the elements are said to be conforming. If 
continuity is violated, the elements are termed 
nonconforming. 

The approximate solution which the finite 
element method provides for the elastic problem 
is determined by making the functional 


E, = 2 U -[ Xuab —f p udB (17) 


D há Bj 


stationary in the class C, of the elastic fields 
piecewise defined by Equation 15. 
UC is the strain energy of the element e 


(Du)! H(Du) dD (18) 
. pe 


If the continuity of the displacements across 
the element boundaries is achieved, the first 
term in F, becomes the strain energy of the body 
and F, becomes the total potential energy. The 
finite element method is then nothing else but a 
particularization of Ritz's method, characterized 
by the piecewise definition of the field. 
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The expression of F, can be modified to take 
into account external forces distributed on the 
element interfaces. F, becomes then 


= Zi (U* — q 


e 
D 


XudD —f pu dB) (19) 
d BI 


B denotes the portion of the boundary of 
element e where prescribed external forces are 
applied. No reactive forces are to be considered. 

In order that the last term of Equation 19 can 


be calculated, the external forces distributed on 
the interfaces must be shared between the pair 


of elements which contact along each interface. 
We assume those forces to be equally distributed 
between the elements in contact.* The magnitude 
of p to be assigned to each element is thus half 
of the distribution density of the force applied 
to the interface. 

Introducing Equation 12 in Equation 19, we 
obtain 


Km | (0 4º)T H (0 fº) dD (21) 
De 
= [ qe” X dD + [ ("| dB 
pe Bi 
(22) 


K“ is the stiffness matrix cf the element. 0º 
is the vector of the external generalized forces 
applied tho the element. 


* This assumption is necessary only if conformity is 
not achieved. 
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Introducting Equation 14, we obtain 
1 T 
E — 2 EH Kº T') 4, — 
e 


(23) 


= 28" ( sr 0º] 


Making 


esp =D [(DF*)” H(D Fº)dD 


(24) 
Te 810" = 5( [FU XaD + 
e [e] = 
De 
+ [ Fº p dB) (25) 
D,e 
there results 
E SF T 
FP, = > q, K, — 14, (26) 


The symmetric matrix K, is called the stiffness 
matrix of the structure. Q, is the vector of the 
external generalized forces applied to the struc- 
ture. 

The stationary conditions for F, are obtained 
by equating to zero the derivztives of F, with 
respect to the mutually independent parameters 
Ji: It results in the system of linear equations 


K, 4, =0, 


n 


(27) 


Equation 27 permits the determination of the 
displacements whenever the structure is stable, 
ie. whenever the prescribed generalized dis- 
placements are enough to prevent rigid body mo- 
tion. Then, indeed, the columns of |[D F“] are 
linearly independent of D and, as H is definite 
positive, K, is nonsingular and also definite posi- 
tive. 


Functional F, can be expressed as 


1 
a 1 Ki —4 K, Jan = 
1 


go: 2 [( “Can y K, (a, Quo) E Van K, Jan | (28) 


in which q,, is the solution of Equation 27. q,, 
corresponds thus to the approximate solution. 
As K, is definite positive, the first term in 
Equation 28 is positive unless q, equals q... 
This proves that the approximate solution mini- 
mizes Fin C,. 
By virtue of Equation 27, 2Uº can be expres- 
e 
sed in terms of the external generalized forces 


applied to the structure. Indeed, 


ai 7 = k (29) 
2 


in which Q is the vector of the external forces 
corresponding to q. 

The same magnitude can also be expressed in 
terms of the external generalized forces applied 
the elements. Indeed, introducing Equation 25 in 
Equation 29, there results 


Su=ia(ST = 
e 


e 


1 e =. 
RAT id =l LO (mi 
e 2 e 
Equation 27 makes it possible to analyse the 
structure by the displacement method. 
The force method may also be used. It starts 


from the reduced equilibrium conditions, which 
are expressed in terms of the vector 


(e = Kº q“ (31) 
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of the total generalized forces acting on each 
element. Multiplying both sides of Equation 31 


by Tº and summing, there results 
é! £6 EM 
A ia Q — > ie Kº q (32) 
e e 


Introducing Equation 14, we obtain 


s TT q — ) Te” Ke Fº qu=K,q, (23) 


and, finally, using Equation 27, 


ST ot=h (34) 


Equation 34 expresses the reduced equilibrium 
conditions. It must not be confused with Equa- 
tion 25. 0º includes indeed not only the gen- 
eralized external forces acting on the element 
(which are the only ones contained in Qº) but also 
the interaction forces between elements and the 
generalized reactions corresponding to the pre- 
scribed displacements. 


SECTION IV 
COMPLETENESS CRITERION 


Completeness of a sequence of families C, with 
respect to a given set C has a meaning provided 
we can compute the distance between any ele- 
ment u, of each family and any element cof C. 
The sequence of families is said to be complete 
with respect to C if it is possible to find, for a 
specified =: >0, an integer N, such that in each 
family with order n >> N there exists an element 


Ucn Which satisfies the inequality 


à (Maj, O) < + (35) 


where c is any element belonging to €. 
The following definition of distance is chosen 
as a basis for our discussion : 
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d(n,, O =1/ 2 f [o mit so | ulo (uº — u ad 
e pe o | 


in which u” and U,. denote, within D*, the dis- 
placement fields corresponding to u, and c. 

The expression under the square root is no- 
thing else than the double of the sum of the 
strain energies corresponding to the difference 
of the two fields u, and c. Such expression is 
never negative and cannot vanish unless the two 
fields coincide (*). It would be easy to show that 
other requirements are also satisfied which make 
Expression 36 a proper definition of distance. 

A general criterion for completeness will be 
stated and justified in this section. Such criterion 
essentially does not differ from the one which 
was presented by Baseley et al. (Reference 5). 

Let C be the set of the compatible elastic 
fields whose displacements have continuous and 
bounded second order derivatives within each 
element. 

We wish to demonstrate that completeness 
with respect to € will be obtained if 

a) the general analytical expression for u; 
within element e is a polynomial with a number 
of arbitrary parameters equal to the number of 
unit modes corresponding to the element, 

b) The terms of degree higher than the first 
can vanish regardless of the values taken by the 
constant term and the coefficients which affect 
the linear terms, 

c) the constant term and the coefficients which 
affect the linear terms are completely arbitrary. 

We remark that, once these conditions are res- 
pected, the displacement component u; or its first 
derivatives cam take up any arbitrary value 
throughout the element, if suitable values are 
ascribed to the parameters. 

The displacements u, of any field c belonging 
to C can be represented inside D“ by the 
following Taylor's expansion 


u=u, (0) +u,, (0) (x,—x)) + 
l 0 o 
+ T Uijx (O)(x,—x,) (x, — x) (37) 


(*) The structure is assumed to be stable, i. e., such 
that the prescribed displacements are enough to prevent 
rigid body motion. 


TECNICA N. 403 


(36) 


O and O, are points in e. O, depends on 
the coordinates of the point where u, isto be 
determined. 

Let us consider now the displacement field 
with components 


Uj=u; (O) + u,,; (O) (x; —x;) (38) 


within D€, which we call tangent field to u at O. 

As all the second-order derivatives are supposed 
to be bounded inside D*, Equations 37 and 38 
yield 


E d 
u, —u; | < ri tra (39) 


in which V, is an upper bound for all the second 


derivatives and é is the maximum diameter of 
element e. d is the total number of second order 
derivatives. 

By considering similar expansions for the first 
derivatives of u,, it is also possible to derive the 
inequality 


Ui,j— Um, | <dV, (40) 


As operator D involves derivatives of the 
first order, we have 


[D(u-u$) | H[D(u-ut) |< vs (Ie)? 


(41) 


for IS sufficiently small. V) is a positive number. 

The distance between the elastic field c, cor- 
responding to U, and u, corresponding to U,, is 
thus 
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d(u, c)= V a [ |D(u—uf) [Em [D(u-ut) | dD < V S Vo (ep — 


“p 


/ 2 pr 
& 1 Vil é Dº =1W'V E p 


in which V, and | denote the maximum values 


of Vo and /º in the whole set of elements. 

Equation 42 means that the distance between 
u, and c tends to zero when n tends to infinity, 
that is, when the size of elements deçreases 
indefinitely. 

Consider now atype of finite element gener- 
ating a sequence of families of fields, C,, whose 
completeness is to be investigated. 

Let u,, be the field of C, whose displacements 
take, at the nodes, the same values as the 
displacement of c. uu, corresponds to u,, 
within D*, 

Suppose the general criterion to be satisfield. 
u, is thus one of the displacement fields which 
can occur within the finite element. Let such field 
correspond to values q,; Of the parameters, i. e., 


U, = f“ q, within Dº (43) 
On the other hand, 
ui = 6º qf (44) 


From Equations 43 and 44, we obtain 


uf-us [=| ff (95, — 98) | vio DS 
(45) 
On the other hand, 
e e e 
Ut, k — Ugi,+| E | a (q: dg) , 
within D* (46) 


But, by virtue of Equation 13, 


8 | TS 
o é 
gem 
Ra 
o 
o) 
no 
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memos 


(42) 


As the absolute dimensions of the element do 
.. . . 4 e 

not appear explicitly in the functions y ;, and 
Xu , these functions remain bounded as the size 
of the element decreases. Assume the moduli of 
all those magnitudes remain below a positive 
number, V.. 

Then 


(48) 


| | V. | 
e e à e e 
Hg Ley | ES — 2 q]; — o; | (49) 
j 


within Dº. 

On the other hand, as u, andu, take the same 
values at the modes, and the parameters q; are 
the nodal values of the displacements, Equation 
39 permits to write 


| e e d ev: 
dy — Ig | Ss Vl p (50) 


Introducing Equation 50 into Equations 48 and 
49, we obtain 


VVNCd o. 
Rr —uç | TT (19 (51) 
e e V, V;Nºd e 
Di Us | (52) 
| ti Cl 21 


Nº being the total number of parameters cor- 
responding to element e. 

Equations 51 and 52 hold even if the first de- 
rivatives of u,, are discontinuous in D”. This is 
an important remark because sometimes (Refer- 
ence 7) the element itself is considered subdi- 
vided into parts and the displacement field admits 


different analytical expressions within each part. 
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The similarity between Equations 51 and 52 
and Equations 39 and 40 allows a jump straight 
to the inequality 


d (ias) Vy, 2 D (53) 


in which u,, represents the piecewise defined 


l 
elastic field with displacements us; within ele- 
ment e. V, is a positive number. 

Combining Equations 42 and 53, we obtain 
finally 


d Mena c) < d (as Ut) a 


(54) 
q q (U in E a VD (Ny v, + Vv,) 


Equation 54 means that the distance between 
U.n and c tends to zero with the size of the 
largest element, so that, as c is an arbitrary ele- 
ment of C, the completeness proof is finally 
achieved. 


SECTION V 
CONVERGENCE DISCUSSION 


Consider any type of finite element which can 
generate a sequence (C, | of families of fields 
complete with respect to C. We wish to investi- 
gate if the sequence of approximate solutions, 
fuan!+ obtained by minimizing F, in each family 
C, converges to the exact solution, which is 
assumed to belong to €. 

We know already that completeness implies 
convergence to the exact solution if it is associa- 
ted with conformity. The finite element techni- 
que becomes then indeed a particularization 
of Ritz's method, in which conformity is a suf- 
ficient condition for convergence. 

Let u,, be the field in C, whose displace- 
ments take at the nodes the same values as the 
displacements of the exact solution. 
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As the exact solution, u,, belongs to C, and 
completeness with respect to C is ensured, it is 
possible to determine N such that, for n > N, 


d (uçiuc)<* (55) 


: being a positive and arbitrarily small number. 
As F, is a continuous functional, we can find £ 
such that 


E (Men) = E, (u5) + & (56) 


* being also positive and arbitrarily small. 
As u., belongs to C,, and u,, minimizes F 
in C,, 


n 


E Ma) SF, (ui) (57) 
and 


F, (un) <F, (u,) + * (58) 


The approximate solution u,, equilibrates a 
system of external forces which generally differs 
from the system of external forces actually 
applied to the body. Let u, denote the compatible 
field which equilibrates the same external forces 
as U,n» i.e., the exact solution corresponding to 
those forces. 

We assume that u, belongs to C, i.e., that the 
second derivatives of the corresponding dis- 
placements are continuous and bounded within 
each element. This assumption will later be 
discussed. 

Let u,, be the field in C, whose displacements 
have the same values at the nodes as the dis- 
placement corresponding tou,. Asu, belongstoC, 
and the sequence (C, | is complete with respect 
to C, it is possible to find N, such that, for 
n>N, 


d (us, Upn) << (59) 


=” being a positive number, arbitrarily small. 
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We can also find N, such that the first deriva- 
tives of the displacements (see Equations 40 
and 52), and thus the difference of the stresses, 
are smaller than «” within each subdomain, that is 


Ra pn | < e (60) 


&” being a column vector with all its elements 
equal to <:”. 

Let us consider now the generalized forces 
connected with the systems of external forces 
respectively corresponding to u, and u,,. 


By virtue of Equation 22 


T 
O; — 0;= [6º 


pº 


(X; — Xº) dD + 


+ (PÉ — pS) dB (61) 
Bj 


in which X; and p, correspond to u,, and X/ 
and p; correspond to u,. 

The body force term tends to be of higher 
order than the second term in Equation 59, as n 
tends to infinity. The difference of the generalized 
forces tends then to be expressed by 


==) (pp — Pp) dB (62) 
e 


By virtue of the continuity of the stresses 
within Dº, Equation 60 is also valid on B". Thus, 


Pp — Pa ES (63) 


and 


e! 9 dB (64) 
Bj 
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On the other hand, by virtue of Equations 18 
30 and 36, the distance between u,, and u,, is 
given by 


2 tê = 7 E = 
(65) 


2 
d(u,, up) | = 


As the external forces which equilibrate u, 


te 
are the same, Q 


“ coincides with 05. Thus, 


andu,, 


e7 
q, E LR 


d (vans Up) <1/ 4 “Sho 


fls 
B| 


Equation 66 shows that the distance between 
an and u,, can be so small as we wish. 
Combining 66 with 59, we conclude that it is 
possible to find N, such, for n > N,, 


u 


SE Edom (67) 


““” being a positive arbitrarily small number. 
As F, is a continuous functional, it is then 


possible, given E , to determine <=“ such that 


F, (u,)=F (us) + + lv (68) 


As u, belongs to C, the theorem of the total 
potential energy permits to write 


E tus) > E tu) (69) 


and thus 


E (u,) <E, (ua) & e" (70) 
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Combining Equations 56, 57 and 68, there 
results 


É Cu SE (ma SEA E É e (71) 
and thus 
E (uc)=E, (ug) + e (72) 
in which 
0o< <p (73) 
But, as u,, and u,, both belong to €,, Equa- 


tion 28 permits us to write 


E (U an) ame (Uan) = (Gen — Gan)” 


Combining Equation 74 and 55 we obtain 
finally 


d (Usa u,) <d (Uan> Uen) E 


+d(umu)<Ce+V2 (75) 


which shows that fu,, | converges to u,. 
It remains to investigate if u, belongs to €. 
Our reasoning will be based on a theorem 
which is known to be valid (Reference 8, p. 349) 
for linear elliptic equations of the form 


2 Ui + É, U,; + yu a É (76) 


Such theorem states that the partial deriva- 
tives of order up to (m+2) of u satisfy a 
Hólder condition with exponent « (0<z<1) 
in every bounded subdomain D“ with closure in 
a closed domain (*) G, whenever all the deriva- 
tives of order up to m of f and the coefficients 


(*) An upper dash denotes a domain and its 
boundary, i e., a closed domain. 


TECNICA Nº 45 


Ko (Men — Us) = [A (Ugo San) = 


aj b, and c satisfy the same Hólde: condition 
in G. 

The displacements will thus be Hólder conti- 
nuous, and therefore have continuous bounded 
second order derivatives in a closed subdomain Dº, 
whenever the body force density is Hólder 
continuous in D*. 

u, is the compatible field which equilibrates 
the same external forces as u,,. Such forces are 
of two kinds: body forces distributed within 
each subdomain, and forces distributed on the 
subdomain interfaces and on B. 

u, will belong to C, that is, the second-order 
derivatives of the corresponding displacements 
will be continuous and bounded within each 
subdomain D”, if the body force density corre- 
sponding to u,, is Hôlder continuous within DS. 

The problem now consists in knowing whether 


Td 


(74) 


the body force density corresponding to u,, is 
continuous and bounded (and therefore Hólder 
continuous) within D*, no matter how large is n. 
This is namely the case if the type of the 
element is such that the body force density 
vanishes or is obliged to a prescribed bounded 
and continuous variation wthin each element, 
no matter the values of the corresponding gen- 
eralized displacements. 

If such a situation arises, there is no doubt 
that the completeness criterion is also a conver- 
gence criterion. 

Our reasoning can be adapted to cases (Ref- 
erence 7) in which the elements are subdivided 
into parts and the allowed fields have different 
analytical expressions within each part. The 
body force density is then generally not conti- 
nuous within the element taken as a whole. 
Convergence will however be proved if each part 
is treated as a separate element. 


SECTION VI 
SIMPLIFIED THEORIES 


The results which have been established can 


be generalized to cover a much broader class of 
problems (Reference 2). The simplified theories 
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which result from neglecting the transverse 
shear deformation (see Section 2) belong to 
such class. 

The convergence criterion respecting the sim- 
plified theories can however be established in 
a more direct way, if the results are used which 
were obtained for the corresponding theories in 
which the transverse shear deformation is not 
neglected. 

It was proved in the last Section that conver- 
gence will be obtained if completeness with 
respect to C is achivied and if the body force 
density components remain continuous and 
bounded within the elements, no matter of 
their size. 

Completeness with respect to C is not however 
a necessary condition for convergence. Com- 
pleteness with respect to a subset CCC 
will indeed be sufficient, if the exact solution 
belongs to C”. 

Assume indeed that all the displacement 
components of every field in C” verify a certain 
set of linear relations 


Ru==0 (77) 


in which R may be a differential operator. 

Consider now polynomial expressions for the 
displacements, within a general finite element, 
and assume that the completeness criterion with 
respect to C is verified. Let C, be the family of 
fields piecewise defined by such polynomial 
expressions made compatible by reduced compa- 
tibility conditions. 

The introduction of Equation 77 implies a linear 
relationship between the coefficients of the 
polynomial expressions which generally involves 
the constant and linear terms. Such terms cease 
thus to be arbitrary and completeness with 
respect to C is therefore destroyed. Completeness 
with respect to C' will however still be main- 
tained if the relations between the coefficients are only 
those which result from Equation 77. Let C, be the 
family of fields piecewise defined by polynomial 
expressions obeving such requirement and made 
compatible by the same reduced compatibility 
conditions as the fields of C. 

The total number of independent coefficients 
corresponding to C, is larger than the total 
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number corresponding to C,. As the number of 
independent cocefficients equals the number of 
nodal displacements, C, corresponds to a larger 
number of nodes and nodal displacements than 
C, and, therefore, to a stiffness matrix K, of 
larger order than the stiffness matrix K, corre- 
sponding to C,. Let q, and q, be the vectors of 
nodal displacements corresponding to C, and C,. 

As the reduced compatibility conditions were 
supposed the same, both for C, and C,, C, is 
a subset of C. The reduced compatibility con- 
ditions are assumed such that the rigid body 
motion of each single element is prevented, so 
that matrices K, and K, are both definite positive. 

To each field in C, corresponds a displace- 
ment vector Q, , as well as a displacement 
vector q, . The elements of vector q, may be 
deduced from the corresponding vector q, by 
using the polynomial expressions of the displace- 
ments to determine the supplementary nodal 
displacements. 

Let q,, be the displacement vector which 
minimizes Fin C,. 

The expression of F, for the fields of C, may 
be given, according to Equation 28, by 


1 
Ea (0 — go) K, (O —Uel= 


 W 
au 2 Gun K, Tan (78) 


As C, is contained in C, , the expression of 
F, for the fields of C| can also be given by 
Equation 78. 

F, will be a minimum in €, if the first term 
in Equation 78 is a minimum. But such first term 
represents half of the square of the distance 
between the field corresponding to q, and the 
field u,, corresponding to Q,,. The field which 
minimizes F, in C, is thus the one in C| whose 
distance to u,, is a minimum. 

Consider the sequence (u | of the approximate 
solutions obtained by minimizing F, in each 
family C,. If the body force density components 
remain bounded within each element, no matter 
the size, the sequence fu | converges to the 
exact solution. 


am 
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But the exact solution was supposed to belong 
to C” and, since sequence (C” | is complete with 
respect to C”, there exists a field u,, CC, whose 
distance to the exact solution tends to zero 
when n tends to infinity. 

The distance between u,, and the exact solu- 
tion also tends to zero. This means that the 
distance between u,, and u,, tends to zero. 

As the distance between u,, and the field 
un Which minimizes F, in C”, is the minimum 
distance between u,, and any field in C.,it 
follows that the distance between u,, and u”,. 


am 
cannot be larger than the distance between u 


an 
and u,,. Therefore, u”,, tends to u,,, and thus 
to the exact solution. 

Completness with respect to a set CCC 
which contains the exact solution, and the con- 
tinuity and boundedness of the body force 
density components, were thus proved to be 
sufficient conditions for convergence. 

It is now easy to establish convergence criteria 
for the simplified theories. The reasoning will 
be exemplified with the theory of plates (Ref- 
erence 9), in which the displacements are the 
normal displacement, U,, and two rotations, 
I, and I,. We assume such magnitudes refer- 
red to a system of cartesian orthogonal coordi- 
nates, x ,, X,, X;. The middle plane of the plate 
corresponds to x,==0, 

The equilibrium equation are 


Nisa +tp;=0 (79) 
M, gs a — Ng; + mp=0 (80) 
and the strain-displacement relations are 
e3=U, +, (81) 
K,g= 13,4 (82) 


In these equations, moments are denoted by 
M, , transverse forces by N,;, transverse shear 
deformations by e,, and curvatures by K,, 
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Greek indices can take up the numerical values 
1 and 2. 

Assume the stress-strain relations to be of 
the form 


N,, = Hp eg; (83) 


Map =Hg,ô K.3 (84) 


The equilibrium equations expressed in terms 
of displacement are then 


H, E (U,,ga + Ia4) + p;=0 (85) 


Ho, py ly, 2 y— Hg, (U,,, + L) + mg=o0 (86) 


if the elastic coeficients H,; and H 
stant all over the domain. 

Assume first that the transverse shear defor- 
mation is not neglected. The results obtained in 
Section V are then directly applicable. This 
means that convergence will be obtained if com- 
pleteness with respect to C is achieved and if 
the body force density components remain con- 
tinuous and bounded within the elements, what- 
ever their size. 

We remember that the body force density com- 
ponents, that is, the elements of vector Y, are 
in this case both the applied force and moment 
distribution densities, p,, m, and m,. Equations 
85 and 86 show that such magnitudes will be 
continuous and bounded if the displacements 
(U,, 1, and 1,) together with their first and 
second derivatives, are continuous and bounded. 

Completeness with respect to C will be obtai- 
ned if the displacements are expressed within 
each element by 


aby5 are con- 


(9) 

| =ag+ayX Hay x +P, (87) 
(2) 

b=ag+tayX Tas x +P, (88) 
(2) 

U;=ay+tanX ag x +P, (89) 


in which po denotes a polynomial whose terms 
are at least of the second degree. The coeffi- 
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cients a, , must be arbitrary and the polynomials 
Pp? can vanish for any value of them. 

Assume now that the transverse shear defor- 
mation is negligible, that is, 


e,;=0 (90) 


This is exactly so if the elastic cocfficients H, s 
are supposed unbounded. 

Introducing Equation 90 into Equation 81 
there results 


L=-—U,,, (91) 


which is a particularization of Equation 77. 
Now let 


(2) p 2 345) 
P; =asx, +agX Xotasxo + Ps (92) 


Introducing Equations 89 and 92 into Equa- 
tion 91, there results 


= (3) 


(3) 


L=-— as— ay xy— 2a5X - Pro (94) 


Comparing Equations 87 and 88 with Equation 
93, it becomes clear that a set of linear relations 
has been introduced between the coefficients 2,.. 

Completeness is still achieved, however, with 
respect to the subset C” of the fields of € with 
vanishing transverse shear deformation (*), pro- 
vided the coefficients a. (i=0,..., 5) are all 
arbitrary. Such completeness, together with the 
continuity and boundedness of the body force 
density components, was shown to be sufficient 
for convergence if the exact solution belongs 
to C”, as it is supposed to. 


(*) As C is the set of the fields whose displace- 
ments are continuous and bounded, together with their 
first and second derivatives, C” is, by virtue of Equation 
91, the set of the fieids in which displacement U , is con- 
tinuous and bounded, together with its derivatives of order 
up to three. 
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As (U,.. + 1.) vanishes and the elastic coef- 
ficients H,. are unbounded, Equations 85 and 
86 cannot be used for the discussion of which 
derivatives must be kept continuous and boun- 
ded in order that m,, m, and p, be also conti- 
nuous and bounded. 

Such difficulty may be removed by remarking 
that the exact solution is in this case the limit 
of a sequence of exact solutions corresponding 
to increasing values of the elastic coefficients Hs... 
Each one of such exact solutions is the limit of 
a sequence of approximate solutions [Un] if the 
completeness criterion with respet to C and the 
condition of boundedness and continuity of the 
displacements and their first and second deriva- 
tives are respected. Such conditions must still 
be valid for the limit case. 

As by virtue of Equation 91, the continuity 
and boundedness of the rotations and their first 
and second derivatives are insured by the conti- 
nuity and boundedness of the first, second and 
third derivatives of U,, it becomes clear that the 
convergence conditions for the simplied theory 
of plates are simply the arbitrariness of the coef- 
ficient «,, (i==0,1,...5) and the boundedness 
and continuity within each element of U, and 
of its first, second and third derivatives. 


SECTION VII 
CONCLUSIONS 


The results of the present paper may be sum- 
marized as follows: 


a) A very simple completeness criterion has 
been given with respect to the set C (to which 
the exact solution is assumed to belong) of the 
compatible elastic fields whose displacement se- 
cond derivatives remain continuous and boun- 
ded within each element, whatever the size. 


b) The completeness criterion becomes a con- 
vergence criterion whenever compatibility is not 
violated across the element boundaries (Ritz 
method). 


c) Completeness with respect to C also en- 
sures convergence in any case, i.e., even if com- 
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patibility is violated, whenever the body force 
density components corresponding to the succes- 
sive approximate solutions remain continuous 
and bounded within the elements as their size 
decreases indefinitely. 


d) Continuity and limitation of the displace- 
ment derivatives of order up to two ensures 
continuity and limitation of the body force 
density components, so that convergence to 
the exact solution will be obtained whenever 
completeness is achieved and such derivatives 
remain continuous and bounded within. the 
elements. 


e) The condition of completeness with respect 
to C is not a necessary condition for convergence. 
Completeness with respect to a subset €' CC 
which contains the exact solution is sufficient. 
This conclusion makes it possible to derive 
convergence criteria for the simplified theories 
which result from neglecting the transverse shear 
deformation. 


f) In what concerns plats, am important 
subset of C is the set C' of the fields of € 
with vanishing transverse shear deformation. 
The completeness criterion with respect to 
C" is simply that the transverse displacement 
U, be given by 


U,=agytagyx tax + ax + 


3 
+ay x Xo+asx+P; (95) 


in which the ccefficients a,, (i=0, ],...5) are 
arbitrary. The polynomial PO can vanish for 
any values of the coefficients, and the rota- 
tions I,, and I, must be deduced by using 
Equation 91. 


2) In the case of the simplified theory of 
plates, suficient conditions for convergence are 
completeness with respect to C”, and the conti- 
nuity and limitation, within each element, 
of the derivatives of U,, of order up to 
three, corresponding to the sucessive approxi- 
mate solutions. 
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h) Completeness with respect to C” is by 
iself, however, a sufficient condition for conver- 
gence whenever compatibility (continuity of U,, 
I, and 1,) between elements is not violated. 


It remains to apply these conclusions to actua] 
cases. 

In what concerns two- and three-dimensional 
elasticity, conformity is easy to achieve. Conver- 
gence will then be obtained if the completeness 
criterion is obeyed. Completeness and the condi- 
tion of the body force density components re- 
maining continuous and bounded within each 
element (this is mainly the case if the body 
forces vanishes, as happens frequently) would 
however ensure convergence even if conformity 
were not respected. 

In what concerns plates (simplified theory), 
conformity has not been easy to achieve. If it is 
not achieved, the completeness criterion (with 
respect to C') will be sufficient for convergence. 
If it is not achieved, convergence is ensured 
only when both the completeness criterion 
and the condition of the derivatives of the 
displacements of order up to three remaining 
continuous and bounded within each element 
are respected. 

As the expressions for U. are polynomials 
at least of the third degree, the third deriva- 
tives of U, generally do not vanish and it 
may happen that they increase beyond all 
limit when the size of the elements decreases 
indefinitely. 

This is why the mesh plays a role in 
convergence. Indeed, some types of mesh 
(like the union jack mesh in Reference 5, 
lead to umbounded third derivatives of U,, 
whereas more regular types of mesh lead to 
bounded third derivatives of U,, and thus to 
convergence. 

Is is possible to check the limitation of the 
third derivatives by inspecting the sequence of 
the appróximate solutions, or by examining the 
behavior of groups of elements (Reference 5). 
Such procedures do not allow prediction of 
convergence before any computations have 
been made, but they make it possible to 
know that the sequence converges to the exact 
solutions, even if such exact solution is not 
available. 
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Levoxin supera a protecção anticorrosiva com hidrazina. 
Considerada até hoje a mais eficaz. 


*Levoxin é um novo tipo de hidrazina 
com capacidade de reacção muito mais 
elevada, que proporciona uma segurança 
de serviço ainda maior nas suas insta- 
lações. 

Levoxin protege da corrosão todo o 
circuito de água-vapor mediante extrema 
aceleração da ligação de oxigénio, al- 
calinização, favorecimento da formação 
de camada protectora de magnetite, e 
graças ao efeito inibidor. 

É muito vasto o campo de aplicação 
deste produto Bayer: em caldeiras de 


toda a espécie, tubagens de vapor e 
de água condensada, aquecimentos a 
distância e por circulação, sistemas de 
refrigeração fechados. E também para a 
conservação, em húmido, de instalações 
que vão estar imobilizadas passageira- 
mente. 

Condições económicas e segurança de 
serviço são os factores decisivos que 
falam a favor de Levoxin. Graças ao 
emprego de Levoxin, reduz-se a ocorrên- 
cia de grandes reparações e evita-se 
dispendiosas falhas de funcionamento. 


timbrado, à: 
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TECNICA XVI 


C. D.U. 535.2: 518.5 


TEORIA DA FONTE MÓVEL (CÁLCULO DA DEPENDÊNCIA 
TEMPERATURA-TEMPO) 


RESUMO 


Utilizando um computador IBM 1620, estuda-se a evo- 
lução da distribuição de temperaturas em função do tempo, 
para uma fonte móvel de calor de velecidade e débito calo- 
rífico constantes, em meios infinito e semi-infinito, resu= 
mindo-se os resultados sob a forma de gráficos. Usam-se 
variáveis adimensionais, 


1. NOMENCLATURA 
1.1 — Valores dimensionais 


T — temperatura 
T, — temperatura inicial 
T, — temperatura de referência 


OQ — débito calorífico da fonte 


V — velocidade da fonte 
« — difusividade térmica 
Y — peso específico 

E card y? ad 22 

t — tempo 


(0, x, y, Z) — referencial associado à fonte 


1.2 — Valores adimensionais 


T—-T, 
= ——— — temperatura adimensional 
I, de Fo 
V A 
n = p= — parametro opera- 
E da tório 
vr : 
p= Eni raio vector adimensional 
x 
V Vy . V 
À me Ma TR da ç a us À y, Z adimen- 
Za Za Za 


sionais, respectivamente 
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SYNOPSIS 


Nondimensional temperature-time dependence for both 
infinite and semi-infinite solids heated by a constant velo- 
city moving heat source is put under graphical form. 
A 1620 IBM computer is used, 


Edf a 
ada = — tempo adimensional 


u, v — constantes com valores possíveis entre 
oe-+ 22 


2. INTRODUÇÃO 


Retoma-se o tratamento do regime transiente 
da fonte móvel num espaço tridimensional, utili- 
zando a expressão de distribuição de temperatu- 
ras na forma que parece mais adaptada às apli- 
cações em meio semi-infinito interessadas neste 
fenómeno [2]: 


3. TRATAMENTO MATEMÁTICO 
3.1 — Objectivo 


Resume-se ao traçado das «isotérmicas» 9/n 
no plano $ == o, para vários q. Para tal, pro- 
cede-se ao corte do espaço por planos Y=u 
e ) ==v, determinando-se º/n nos pontos inter- 
secção dos conjuntos de três planos. A gama 
de valores de ) considerada será — 100 << À < 2. 
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32 — Processamento 


Fez-se através de um computador IBM 1620. 
O programa, em linguagem FORTRAN II, repro- 
duz-se na gravura n.º 1. As funções erf foram 
substituídas por um desenvolvimento aproximado 
escolhido entre os que constam de [4] e que, 
envolvendo constantes que «cabem» no compu- 
tador IBM 1620, conduz a resultados com erro 
de valor perfeitamente satisfatório e, principal- 
mente, garante a convergência da segunda par- 
cela da chaveta de (1). Não se transcreve por 
ser de imediata identificação através do programa. 


DIMENSION SIGHA(S0), CANB(GO),RO(60), TETANÇ50) ,PSI(10) 
| REsteio rd 
2 FORMAT(313) 
L FORMAT(F 9,3) 
FORMAT(F 7.3) 
FORMAT(FS.1) 
1 FORMAT(2ZOX bHLAHNB, 9X, SHTETAN,/) 
FURJAT( SHP PY 


| 
IR LHPSl= F5,2) 
12 FORMAT(6HS | 
70 FORMAT(1I7X,F7.3,5X,ElL,o) 
READ 1,A1,A2,A3,P 
READ Z,N,I!,K 
Do 3 |=1,M 
3 READ L,CAMB(1) 
DO 5 J=1,M 
5 READ 6,PSI(J) 
DO 7 L=1,kK 
7 READ O,SIGMA(L) 
BEGIN TRACE 
DO 60 L=1,K 
PRINT 12,SIGMA(L) 
DO 60 J=1,M 
PRINT 11,PSI(J) 
DO 9 1=1,N 
RO(I)=SORT(CAMB( 1 )==24PSI(J)+=2) 
OCO O IS do as SALTA co TEMA L) 
X2=RO(I)/(2.*SIGMA(L))+SIGMAÇL) 
IF (X1)30,30,40 
30 X1=X1=(=1.) 
T=1./(1.+P*X1) 
PI=A1ZT4.A2AT=*24A 3473 
Q=EXP(X|=*2) 
R=P1/Q 
EFCX1=2.-R 
GO TO 50 
CONT INUÉ 
T=1./(1.+P=X1) 
PISAJST-AZOT *24A307 +3 
O=EXP(X10*2) 
R=P1/Q 
EFCXI=R 
50 T2=1./(1.-P*X2) 
PI=ANTZAZETIOAZLA TETO 
Q2=EXP(X2:*2) 
R2=P2/Q2 
EFCX2=R2 
TETA=(EFCXI /EXP(RO(I))+EFCX2=EXP(RO(I 2. "RO(I 
9 TETANLTO TOR RPI AR CRO(I9))/6 RNA 
PRINT 10 
DO 60 |=1,N 
60 PRINT 7O,CAMB(I),TETAN(1) 
END TRACE 
STOP 
END 


ad 


Fig. 1 


«A instrução «begin-end» «frace», incluída no programa, 
será ignorada pelo computador se a chave operatória 
n.º 4 se encontrar na posição «off», Ela permitirá, se 
desejado, acompanhar a evolução das variáveis 
intermediárias». 
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3.3 — Resultados 


Reproduzem-se em gráficos “/n == F (4,1) nas 
gravuras n.0S 2 a 8. Dada a simetria do fenó- 
meno, omitiu-se, por desnecessária, a metade 
inferior dos gráficos. 


4. COMENTÁRIO 


4.1 — Os resultados calculados, para valores 
de q elevados (figs. 2 a 5), conduzem a gráficos 


Fig. 2 


inteiramente concordantes com os inclusos em 
[3] (regime estacionário do fenómeno em variá- 
veis adimensionais). 


2.0) 


e 


19 0, 


-1 -28 


| 


Fig. 3 


4,2— O gráfico n.º 6, para 7 = 5, demons- 
tra que, para menores valores deste parâmetro, 
se verifica a «contracção», na cauda da fonte, das 
«isotérmicas» de menor valor. Neste caso apre- 
senta-se apenas uma figura porque as isotérmi- 
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a | 
5 
E 5 ) 
-06 -04 -o2 02 04 


cas superiores a 0,02 serão coincidentes com as 
do regime estacionário. A tentativa da função 
erfc da primeira parcela da chaveta de (1) para, 


C=50 f 
=10 Y | P 40] 
- 0,01 1 40 
a A qt ce ço 
- 100 -S0 -60 -40 -20 | 0º. . 
Fig. 5 


através do seu argumento negativo, conservar o 
campo de temperaturas em valores vizinhos dos 


estacionários é vencida pelo crescimento inevi- 
tável de ; quand 2 atinge valores de ordem — 35 
a— 50. 


4.3 — As figuras para 7 ==1 mostram agora, na 
cauda da fonte, um compacto grupo de isotér- 
micas, fechando-se todas para À > — 5. Como é 
evidente, as modificações na frente da fonte são 
muito menos sensíveis. As «isotérmicas» supe- 
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Fig. 7 


riores a 5 serão, práticamente, idênticas às esta- 
cionárias. 


5. CONCLUSÕES 


1) Em processos envolvendo fontes móveis que 
interessem o meio infinito, as «isotérmicas» (/n 
serão ainda as mesmas, se forem usadas as 
expressões definidas na nomenclatura. Apenas ao 
separar O de n, para obter os seus valores cor- 
rectos, será necessário dividir ambos por dois. 


| 
A “4 
0=1 y | 
que 
AS 
oó ê 
o 
NL o 
ei -2 | 2 


Fig. 8 


2) Só o conhecimento dos parâmetros parti- 
culares dum processo, bem como da natureza do 
mesmo, poderá permitir definir, de acordo com 
o fenómeno específico a observar, o interesse da 
consideração do regime transiente — pela sua 
duração — nesse processo, 


3) O programa usado permitirá a investigação 
de qualquer 7 particular desejado. 
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UTILIZAÇÃO DE COORDENADAS PARABÓLICAS 
NA ANÁLISE DE TENSÕES DE UM CORPO 
DE FORMA PARABÓLICA 


RESUMO 


Pretende-se com este trabalho deduzir a expressão geral 
do equilíbrio em coordenadas parabólicas, partindo da 
expressão geral do equilébrio em coordenadas curvilineas, 

Justifica-se em linhas gerais o emprego das coordenadas 
purabólicas mo estudo de um corpo crjo contorno tenha 
a forma parabólica, 


NOTAÇÃO 


h — factor de escala 
p — distância 
r—raio de curvatura 
x, Y, Z, — coordenadas cartesianas 
u, v, Z,— coordenadas parabólicas 
(x, *) — coordenadas curvilíneas 
o — tensão normal 


7 — tensão tangencial 


1— INTRODUÇÃO 


Ao fazer-se o estudo da distribuição de tensões 
num corpo contínuo qualquer, um dos primeiros 
pontos que se deve ter em atenção é a escolha 
de um sistema de coordenadas conveniente. 

Um dos factores que condicionam a escolha 
do sistema de coordenadas é a forma que o corpo 
apresenta, uma vez que a simplicidade do estudo 
a realizar está intimamente relacionada com o 
binário sistema de coordenadas — forma do 
corpo. 

Vamos, neste artigo, justificar o emprego de 
coordenadas parabólicas no estudo de um corpo 
com forma parabólica, e deduzir as equações de 
equilíbrio naquele sistema de coordenadas. 
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SYNOPSIS 


À general expression vf equilibrium in parabolic coordi- 
nates is obtained, starting from the general equation of 
equilibrium in curvilinear cooroinates, 

In general terms it is justificd the use of paraboliecadrdi- 
nates in the stress analysis fr a body of parabolic-furm, 


lg os 2olsl 
2 — COORDENADAS PARABÓLICAS |... 


2.1 — Utilização saclas O 


as log eobsb 


Um dos mais importantes condicionamentos 
na escolha das coordenadas a utilizar no estudo 
de qualquer corpo é a forma que o corpo apre- 
senta, sobretudo pela necessidade de definir as 
condições de fronteira ou condições aos limites, 
isto é, os valores que os parâmetros em jogo 


tomam em pontos definidos do contorao do,gorpo. 
Interessa-nos pois trabalhar cem. um sistema 
de coordenadas em que os valores das. variáveis 


tomados na fronteira do corpo sejam aí valores 
constantes. Ora isto consegue-se para, q gaso de 
um corpo com forma parabólica, utilizando um 
sistema de coordenadas parabólicas. : 
As vantagens deste procedimento” são análogas 
às que se obtêm quando se empregam coofde- 
nadas cartesianas no estudo de um corpo para- 


lelipipédico ou coordenadas esféricas. no estudo 


, 
de um corpo esférico. sv 20 as8luoÍso sup 


| JET1, 6 OED,.PON SUP 
22 — Dedução das equações dé equilibrio“. 


2.2.1 — Descrição do Sistema de Coordenadag7i222 


O sistema de coordenadas, considerado é com- 
posto por duas parábolas co-focais com os mes- 
mos eixos, conforme se ilustrarna''fig. +. vb 


mm. te TO -m 
«A família de curvas de v=-c eu=c são 
cilindros parabólicos, cujos traços no plano xy 
se intersectam perpendicularmente e que têm os 


Fig. 1 


vértices no eixo dos xx mas em lados opostos 
. 4 * te — 

relativamente à origem. z==c são planos, para- 

lelos ao plano xy. 


2.2.2 — Dedução matemática 


Os valores de x e y em função de u e v são 
dados pelas expressões seguintes: 


| 
= 
< 


y 


A equação de equilíbrio em coordenadas curvi- 
líneas para a direcção « é dada pela expressão (1) 
e torna-se necessário agora estabelecer as relações 
entre dz, dÉ, r, e r; em termos do novo sis- 
tema de coordenadas. 


da + dB + r; a ii (1) 


Para determinar os valores de d x e d É teremos 
que calcular os valores dos factores de escala 
que nos dão a transformação de da e d2 em 
termos de du e dv. 

Assim : 


dx = udu — vdv 


dy = udv + vdu 


130 


e (dp) = (dx) + (dy) 
= (u? + v) (du) + (u” + vº) (dv) 


logo os factores de escala são: 


h, ni Vu + vi 


ig = A 4 É 


De modo a obter os valores de r, er :, temos 
que determinar os seus valores através da expres- 
são geral do raio de curvatura dado por 


d” y 
E d x? 
= , “2 
dy 

1 hat 

[+ (8) 

ESA dy 4 dy dy 
dy dy du | du — du du 
dx du dx dx dx 4 Ox dv 

du du dv du 

dv 

vtu— 

É 

dx div 

E sy 
du 
Para v == const. > MP. 
dx u 


Fazendo agora 


= D=t(u,v, E 
dx du / 
dº d 
= (8) 
dx dx 
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O ad di 
dx dv 
u—v— 
du 
dv dv dv > dv v dv 
o (du FR + u pe du + dv) (uv o) (ua +) (du-dy ny E du ) 
(u- vao) (udu — vdv) 
' du 
Fazendo 
- dv 
V' = — 
du 
d'y — Cuví + 2vº) (u—vv') — (uv +v) [1— (v)*— vv” ] 
dx” (u-v é) 
du 
= 


v—uvvv +auv—zv(v) uv —vtu(v) + v(v)tuvv vo + viv” 
(u—vv) 
uv” +uv—v(v)-v+u(v) +v?v” 


(u—vv) 


Para v == const. 


Deste modo a expressão final será: 
3 
EE is 1º day 1º 
dx u” Vu? +v' du Vu +v” dv 
e 
V Suu — Soy > Cuv 
RO + ii 32 = 0 
1 == u o ; 
= v 
= x a (1 E ) (a + A ) 
(1+5) 
u u v 
ás 3 . 3 
v u 
52 
v? 
1 q? Para a direcção £ 
= 
v 
a 1 54 1 Tou 
Vu +v? dv Vu? + v” du 
512 
uí v 
(a EE x ) + TO UR = (3,9 — Tuy) — 
Sã La t-— 
3 2 
v v 
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Simplificando : 


Direcção “u”: 


16b 


dO uu 


du 


dTuv 


dv uty 


ut v? uv 


2v 


u 


(Tuu” A" NR 


= O 


(2) 


Direcção “v”: 


dTw dTeu 
-- ) ) ; Cove 9 
dv du 3 ( vv dt) 
do (3) 


= 2 v 
uv u 


Estas expressões (2) e (3) são equações de equi- 
líbrio para as direcções “u” e “v” num sistema 
de coordenadas parabólicas e deduzidas a partir 
das equações de equilíbrio em coordenadas curvi- 
líneas na sua forma mais geral. 
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A aplicação de depósitos de cromo no campo industrial, pelo sistema de electrodeposição, é 
relativamente recente. Inicialmente estes tratamentos eram sobretudo utilizados com fins decorativos, 
depositando uma ligeira camada de crómio sobre um revestimento geralmente de níquel, para se obter 


um aspecto brilhante e resistente à corrosão. Tais depósitos são macios e de pouco valor para a uti- 
lização industrial. 


As propriedades do cromo, quando depositado em condições perfeitamente controladas, são tais 


que levaram à sua larga aplicação na moderna tecnologia mecânica. É o processo conhecido pelo nome 
de CROMAGEM DURA. 


As propriedades que oferece o CROMO DURO, tanto no tratamento de peças novas como 
nas recuperadas, são principalmente as seguintes: 


— elevada dureza (até cerca de 1000 Vickers). 
— baixo coeficiente de atrito. 

— alta resistência à corrosão e ao desgaste. 
— ausência de distorção. 

— possibilidade de aplicações locais. 

— aderência perfeita ao metal de base. 


Existem três tipos de cromagem dura — ligeira, média e espessa. Para cada aplicação deve esco- 
lher-se o tipo mais conveniente. 


Em qualquer dos casos, e mesmo que as peças sofram rectificação posterior, é necessário que 
sejam rectificadas e poldas antes de se efectuar o revestimento. 


O nosso Departamento Técnico está preparado para prestar esclarecimentos, estudar e propor 
as melhores soluções para os problemas que lhe forem submetidos. 
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CONCEITOS E PROPRIEDADES TOPOLÓGICAS DO 
MÉTODO DE CÁLCULO DE CIRCUITOS ELÉCTRICOS 
PELA TEORIA DOS GRAFOS 


RESUMO 


Baseado nos conceitos da teoria dos grafos estudam-se 
as propriedades fundamentais relacionadas com à estru- 
tura, planaridade, singularidade, corte e dualidade dos 
grafos. Para o cúlculo de circuitos eléctricos definem-se as 
matrizes de incidência que interessam nos diferentes méts- 
dos: dos nós, das malhas, dos cortes e dos anéis, 


1. INTRODUÇÃO 


A topologia (do grego «lugar» e «ciência») é 
a parte da geometria que, abstraindo das gran- 
dezas mensuráveis de natureza física ligadas à 
estrutura geométrica, tais como ângulos, compri- 
mentos, secções, massa, resistividade, permeabi- 
lidade, etc., se ocupa do modo de ligação entre 
os diferentes elementos quaisquer que sejam as 
deformações contínuas sofridas por esses ele- 
mentos. 

A esfera e o cubo, por exemplo, são topolo- 
gicamente iguais visto que em cada figura o 
interior é separado do exterior por superfícies 
topológicamente de igual valor, pois em ambos 
os casos se trata de superfícies simplesmente 
conexas (que se podem reduzir a um ponto por 
deformação contínua). 

Qualquer figura, por exemplo uma circunferên- 
cia, traçada sobre uma membrana de borracha é 
topológicamente a mesma que toda a figura 
resultante de qualquer deformação da membrana 
(sem a quebrar). Um rectângulo cujos lados e 
diagonais sejam fios de borracha é igual sob o 
ponto de vista topológico a qualquer figura com 
os mesmos nós e igual número de ramos mas 
com os lados deformados. 

O caminho através de um labirinto ou o 
entrançado dos fios de renda ou ainda os nós dos 
cabos em marinharia são questões topológicas. 

A electrotecnia, no cálculo de complexos cir- 
cuitos eléctricos, serve-se da parte da topologia 
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SYNOPSIS 


The basic conceptions of graph theory are presented and 
the following properties are stwlied: atructure, planarity, 
singularity, cut and duality of graphs. For electrical 
network caleulation the node, mesh, cut and loop incidence 
matrices are defined and discussed, 


que estuda as configurações lineares, ou seja da 
teoria dos grafos, para obter um tratamento mate- 
mático simples, através de um método exacto e 
seguro, dando resposta às seguintes questões 
fundamentais: Quantas equações são necessárias e sufi- 
cientes para calcular um dado circuito eléctrico? Como se 
consegue definir um sistema de equações independentes ? 

Uma vez obtido o sistema de equações a sua 
resolução processa-se pelos métodos da álgebra 
linear. Neste sentido o cálculo matricial permite a 
resolução de sistemas complexos de um modo 
elegante e simplificado. 

Para se compreender contudo a aplicação do 
método de cálculo introduzem-se as noções 
fundamentais da teoria, repartidas em definições 
(que se indicam abreviadamente por DEF.) 
e respectivas propriedades (referenciadas pela 
abreviatura PROP.), de um modo sistemático 
esclarecendo-se gradualmente os conceitos essen- 
ciais. 

O presente trabalho foi elaborado em 1969 
a propósito do estudo das redes eléctricas na 
Cadeira «Aplicações de Electricidade II» do curso 
de Engenharia Electrotécnica da Universidade 
de Luanda. 


2. CONCEITOS DA TEORIA DOS GRAFOS 
2.1 — Estrutura dos grafos 


DEF. 1 — Grafo é toda a figura formada por nós 
e ramos (fig. 1-a). 
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DEF. 2— Nó de um grafo é qualquer ponto de 
ligação de pelo menos dois ramos (pontos 1 a 10 da 
figura 1-a). 

DEF. 3 — Nó principal é um nó que liça mais de 
dois ramos (nós 2 a 10 da fig. 1-a). 

DEF. 4 — Nó secundário é um nó que liga apenas 
dois ramos (nó 1 da fig. 1-a). 

DEF. 5 — Nó origem é um nó arbitrário tomado como 
referência. 

DEF. 6 — Par de nós é o conjunto de dois nós 
quaisquer, 

PROP. 1 — Qualquer nó de um grafo constitui 
um par de nós com o nó origem. 

DEF. 7 — Ramo de um grafo é o elemento que 
liga dois nós. 

PROP. 2 — Eliminando os nós secundários o 
númerototal de ramos reduz-se de igual quantidade. 

PROP. 3 — As extremidades dos ramos são 
nós, isto é, as extremidades de um ramo não 
fazem parte desse ramo (um ramo é aberto). 

DEF. s — Sub-grafo de um dado grafo é um 
erafo particular que contém apenas parte dos nós e ramos 
do grafo total (fig. 1-b). 

DEF. 9 — Sub-grafo complementar a um sub- 
-grafo de um dado grafo é o restante sub-grafo no 
grafo total (fig. 1-b). 

DEF. 10 — Grafo ou sub-grafo linear é um grafo 
que contém apenas nós secundários excepto os nós extre- 
midade (fig. 1-c). 

PROP. 4 — Exceptuando os nós extrémidade 
de um grafo linear não fechado todos os restan- 
tes nós ligam apenas dois ramos. 

DEF. 11 — Percurso entre dois nós de um grafo 
e um sub-grafo linear tendo esses nós como extremidades 
(fig. 1-d). 

DEF. 12 — Anel de um grafo é qualquer sub-grafo 
linear fechado do grafo (fig. 1-e). 

PROP. 5 — Num anel todos os nós ligam dois 
ramos. 

PROP. 6 — Retirando um ramo a um anel 
obtém-se um percurso. 

DEF. 13 — Árvore (!) de um grafo é um sub-grafo 
que contém todos os nós do grafo ligados por ramos sem 
a formação de anéis (fig. 1-f). 

PROP. 7 — Num grafo existe pelo menos uma 
árvore. 

PROP. 8 — Cada ramo de um grafo pode fazer 
parte de uma árvore. 


(') Esta noção é por vezes referida por árvore com- 
pleta |1). 
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PROP. 9 — Uma árvore contém em geral nós 
principais e secundários. 

PROP. 10 — Numa árvóre existe apenas um 
percurso entre cada par de nós. 

PROP. 11 — Não existe nenhum anel numa 
árvore. 

DEF. 14 — Co-árvore de um grafo é o sub-grafo 
complementar de uma árvore do grafo (fig. 1-f). 

PROP. 12 — Um grafo é definido por uma 
árvore e respectiva co-árvore. 

DEF. 15 — Haste de uma árvore é qualquer ramo 
de uma árvore do grafo. 

PROP. 13— Quando se suprime uma haste 
a uma árvore esta divide-se em duas árvores 
(considerando um nó isolado como uma 
árvore). 

DEF. 16 — Árvore linear de uma árvore é um 
percurso entre dois nós extremidade contendo todos os 
nós principais da árvore (fig. 1-g). 

DEF. 17 — Ramificação de uma árvore é uma 
haste que liga um nó principal a um nó extremidade da 
árvore excluindo a correspondente árvore linear (fig 1-g). 

PROP. 14 — Uma árvore que seja um grafo 
linear não tem nenhuma ramificação (pois não 
possui nenhum nó principal). 

DEF. 18 — Corda de um grafo em relação a uma 
arvore é qualquer ramo da respectiva co-árvore. 

PROP. 15 — Qualquer corda não pertence à 
árvore do grafo. 

PROP. 16 — Num grafo existe sempre igual 
número de cordas independentemente da árvore 
escolhida. 

PROP. 17 — Um grafo pode definir-se por uma 
árvore e pelas relativas cordas. 

DEF. 19 — Anel-base de um grafo é um anel que 
contém apenas uma corda (fig. 1-h). 

PROP 18 — Num grafo existem tantos anéis- 
-base quantas as cordas. 

PROP. 19 — Juntando uma corda a uma árvore 
obtém-se um anel-base, que é a malha contituida 
pela corda e pelo percurso única na árvore entre 
as extremidades relativas à corda. 

DEF. 20 — Grau de conexidade c de um grafo 
é o número de anéis-base (ou de cordas) do grafo. 

PROP. 20— O grau de conexidade é indepen- 
dente da árvore escolhida, sendo únicamente 
função do número de ramos r e do número de 
nós n do grafo. 

PROP. 21 — Uma árvore é um grafo simples- 
mente conexo (uma vez que não contém nenhuma 
corda). 
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PROP. 22 — Juntando uma corda a uma árvore 
forma-se um grafo de conexidade unitária (por 
se formar um anel-base). 


f) 9) h) 


Fig. 1 — Elementos estruturais dos grafos. 
a) Grafo com n=-I0 nós (1 secundário) e 
r = 16 ramos; 
6) Sub-grafos complementares ; 
c) Grafo linear; 
dy Percurso; 
e) Anel; 
f) Árvore e co-árvore; 
£) Árvore linear e 4 ramificações; 
h) Anel-base, 


2.2 — Planaridade 


DEF. 21 — Transfiguração de um grafo consiste 
na alteração da sua forma por deformação contínua sem 
perda de conexidade (fig. 2-a). 

DEF. 22 — Grafos semelhantes são os que por 
transfiguração podem assumir a mesma configuração 
geométrica (fig. 2-a). 

PROP. 23 — É condição necessária mas não 
suficiente para que dois grafos sejam semelhantes 
que tenham igual grau de conexidade. 

DEF. 23 — Cruzamento de ramos de um grafo 
é um ponto onde se cruzam dois ou mais ramos sem que 
seja um nó, ou seja, é um ponto onde se verifica simples 
sobreposição dos ramos (fig. 2-b). 

PROP. 24 — O cruzamento de dois ramos quais- 
quer pode evitar-se pela transfiguração de um 
dos ramos (fig. 2-b). 
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DEF. 24 — Cruzamento autêntico ou encruzi- 
lhada é um cruzamento que não se pode evitar pois a 
transfiguração conduz ao deslocamento do cruzamento 
para outros ramos (fig. 2-c). 

PROP. 25 — Uma encruzilhada não é ligada a 
dois ramos determinados mas pode deslocar-se 
para diferentes ramos por transfiguração do 
grafo. 

DEF. 25 — Graáfo planar é um grafo que se pode 
representar num plano ou numa esfera sem cruzamentos 
(fig. 2-b). 

PROP. 26 — Uma árvore é um plano planar. 

DEF. 26 — Malha de um grafo planar é todo o 
ânel que numa representação sobre uma esfera não contém 
ramos no interior, 

DEF. 27 — Malha Interior de um grafo planar 
é um anel que numa representação no plano não contém 
ramos no interior (anel 1-2-3-1 no segundo grafo 
da fig. 2-b). 

DEF. 28 — Malha exterior ou cercadura de 
um grafo planar é o anel que numa representação no 
plano não contém ramos no exterior (anel 1-2-4-1 no 
segundo grafo da fig. 2-b). 

PROP. 27 — Pode definir-se malha de um grafo 
planar como todo o anel que, numa representação 
gráfica no plano do grafo, não possui ramos no 
interior ou no exterior. 

PROP. 28 — Na representação de um grafo 
sobre uma esfera nenhuma malha se individualiza 
como malha exterior. 

DEF. 29 — Grafo não planar é um grafo em que 
existem cruzamentos autênticos quando representado sobre 
um plano ou uma esfera (fig. 2-c). 

PROP. 29 — Um grafo não planar contém pelo 
menos quatro malhas (visto que um percurso 
situado inteiramente fora do plano forma pelo 
menos uma malha com a porção do grafo 
situado no plano). 

PROP. 30 — Existe sempre uma superfície 
espacial de conexidade suficiente sobre a qual 
um dado grafo não planar pode ser representado 
sem cruzamentos. 

PROP. 31 — Existe uma superfície espacial de 
conexidade mínima sobre a qual um dado grafo 
pode ser representado sem cruzamentos. 

DEF. 30 — Gênero q de um grafo é o grau 
mínimo de conexidade da superfície orientável sobre a 
qual o grafo pode ser representado sem cruzamentos. 

PROP. 32 — Um grafo planar é de género nulo 
(dado que se pode representar sem cruzamentos 
sobre uma esfera, a qual possui género nulo). 
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PROP. 33 — Um grafo representável sem cruza- 
mentos sobre uma superfície toroidal (que tem 
género g==1) é um grafo de género unitário. 

DEF. 31 — Grafos elementares não planares 
são o pentágono completo (definido pelos cinco lados 
e pelas cinco diagonais) e o hexágono completo 
(definido pelos seis lados e pelas seis diagonais). 

PROP. 34— O pentágono completo e o hexá- 
gono completo têm género unitário, pois se podem 
representar sem cruzamentos sobre um toro e 
possuem encruzilhadas quando representados 
num plano ou numa esfera (fig. 2-d). 

PROP. 35 — Teorema de Kuratowski: É condição 
necessária e suficiente para que um grafo seja planar 
que não contenha como sub-grafos nem o pentágono 
completo nem o hexágono completo. 

PROP. 36 — A não planaridade de um grafo 
revela-se quando se atinge um dos grafos elemen- 
tares não planares (consequência do teorema de 
Kuratowski). 


DEF. 34 — Nó singular ou nô de separação 
num grafo conexo é um nó que constitui o único ele- 
mento a um dado sub-grafo e no seu complemento 
(fig. 3-c). 

PROP. 38 — À condição necessária e suficiente 
para que um dado nó seja nó de separação 
é que existam dois outros nós do grafo tais que 
qualquer percurso que os una passe necessária- 
mente por esse nó. 

PROP. 39 — Cada nó interno de uma árvore 
é um nó de separação. 

DEF. 35 — Desdobramento de um nó de sepa- 
ração consiste na reprodução deste nó em tantos nós 
quantos os sub-grafos isolados que pode provocar 
(fig. 3-e.) 

PROP. 40 — Os ramos inicialmente incidentes 
num nó de separação incidem, após o desdobra- 
mento, em dois ou mais nós (tantos quantos os 
sub-grafos isolados do grafo conexo inicial resul- 


tantes do desdobramento). 


| 
HH 


d) : 


Fig. 2 — Planaridade dos grafos. 


a) Transfiguração; 
b) Grafo planar; 


c) Grafos elementares não planares ; 
dj Grafos de género unitário (pentágono e hexágono completos). 


2,3 — Singularidade 


DEF. 32 — Grafo conexo é um grafo em que se 
pode chegar a um nó partindo de qualquer outro nó por 
meio de um certo número de ramos (fig. 3-a). 

PROP. 37 — Entre cada par de nós de um 
grafo conexo existe pelo menos um percurso. 

DEF. 33 — Grafo não conexo é um grafo composto 
de pelo menos duas partes completas separadas (fig. 3-b). 
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DEF. 36 — Haste singular ou haste de separa- 
ção é o ramo de árvore cujas extremidades são nós de 
separação (fig. 3-e). 

PROP. 41 — Desdobrando todos os kh ramos 
de uma árvore esta decompõe-se em nl 
ramos isolados, onde n, é o número de nós da 
árvore. 

DEF. 37 — Grafo separável é aquele que contém 
pelo menos um nó singular, 
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2.4 — Corte 


DEF. 38 - Corte de um grafo consiste na extrac- 
ção simultânea de um conjunto de ramos que provoque 
a sua cisão em dois sub-grafos isolados (fig. 3-f). 

DEF. 39 — Feixe de corte é o conjunto de ramos 
do grafo resultantes do corte. 

PROP. 42 — O corte indica-se graficamente por 
uma linha de separação fechada e os ramos por 
esta intersectados constituem o feixe de corte. 

PROP. 43 — Um feixe de corte contém pelo 
menos uma haste da árvore. 

DEF. 40 — Corte-base é um corte em que se extrai 
apenas uma haste da árvore do grafo (fig. 3-£f). 


Y S 
SZ 
7 as 


b) Cc) 


e os seus ramos (unindo os nós duais das malhas 
contíguas do grafo dado). 

PROP. 48 — O grafo dual de um grafo planar 
é também planar (consequência das malhas do 
grafo dual se poderem representar sobre um plano 
ou esfera). 

PROP. 49 — É condição necessária e suficiente 
para que um grafo admita dual que seja planar 
(pelo que se pode definir um grafo planar como 
todo aquele que tem um grafo dual). 

PROP. 50 —-Se um grafo G' é dual de um 
grafo G, então G é dual de G”. 

PROP. 51 — O grafo dual de um dípolo é um 
dípolo (para se obter o dípolo dual pode conce- 
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Fig. 3 — Singularidade e corte de grafos. 


a) Grafo conexo; 


b) Grafo não conexo; 

c) Nó de separação (singular) ; 
d,y Desdobramento ; 

e) Haste de separação (singular). 


DEF. 41 — Feixe de corte-base é um feixe de corte 
que contém somente uma haste da árvore. 

PROP. 44 — Para uma dada árvore de um grafo 
cada haste define univocamente um feixe de 
corte-base. 

PROP. 45 — Um grafo possui tantos feixes 
de corte-base quantas as hastes (ou seja, h== 
= n,— 1). 


2.5 — Dualidade 


DEF. 42 — Dual de um grafo é um grafo em que 
cada nó corresponde « uma malha do grafo dado, cada 
ramo a um ramo do grafo dado e cada malha a um 
nó do grafo dado (fig. 4-a). 

PROP. 46 — Na malha exterior existe um nó 
“dual (exterior ao grafo dado). 

PROP. 47 — Na construção de um grafo dual 
marcam-se os seus nós (um em cada malha inte- 
rior e na malha exterior do grafo considerado), 
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ber-se o dípolo fechado por uma linha; figura 
4-b). 

PROP. 52 — A dualidade topológica caracteriza-se 
pelas correspondências : 


nó «——, malha 
haste «——», corda 
árvore «——, co-árvore 
feixe de corte «——» anel 
ligação em série <-—-., ligação em paralelo (fig. 4-b) 
ligação em estrela «——., ligação em triângulo (fig. 4 a) 


DEF. 43 — Grafo planar orientado é um grafo 
em que cada ramo possui uma flecha indicadora de um 
sentido positivo relativo a uma determinada grandeza 
física (por exemplo, tensão, corrente, energia, etc.) 
(fig. 4-e). 

PROP. 53 — Cada ramo orientado é conver- 
gente para um seu nó extremidade e divergente 
do outro. 
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PROP. 54 — A todo o grafo planar orientado 
corresponde univocamente um grafo planar dual 
orientado. 

PROP. 55 — Grafos duais orientados têm os 
ramos orientados em correspondência com duas 
grandezas físicas duais. 

DEF. 44 — Orientação dual de dois grafos 


duais é definida de modo a que o ângulo entre os ramos 


a) 


PROP. 58 — Numa árvore o número de hastes 
é inferior de uma unidade ao número de nós 
(em virtude da definição de árvore) 


h=n-1 (3) 


e, inversamente, se um grafo tem n nós e n-1 
ramos então é uma árvore (em particular, esta 
relação é válida para um percurso). 


Fig. 4 — Dualidade de grafos. 


a) Grafos duais; 


b) Dualidade de ligação em série e em paralelo; 
c) Orientação dual de malha interior; 

d) Orientação dual de malha exterior ; 

e) Grafos duais orientados. 


duais seja inferior (ou superior, conforme a con- 
venção) a = radianos medindo sempre os ângulos a 
partir dos ramos de um dos grafos para os correspon- 
dentes do outro no mesmo sentido de rotação (conven- 
cionalmente do ramo dual para o ramo inicial no 
sentido directo) (fig. 4-c, d, e). 


3. RELAÇÕES TOPOLÓGICAS FUNDAMEN- 
TAIS 


DEF. 45 — Relação topológica consiste numa 
expressão algébrica entre os elementos de um ou mais 
grafos. 

PROP. 56 — Para se ligarem n nós são precisos 
pelo menos 1-1 ramos 


r> nl (1) 


PROP. 57 — Numa árvore de um grafo o 
número de nós n, é igual ao número total de 
nós do grafo n (em consequência da definição de 
árvore) 


n.=n (2) 
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PROP. 59 — Numa malha e num anel o número 
de nós n, é igual ao número de ramos r, 


=, (4) 
PROP. 60 — Numa árvore o número de nós 


extremidade rn, é superior em duas unidades ao 
número de ramificações R 


ne=R+2 (5) 


pois uma árvore não ramificada possui duas 
extremidades e cada nova ramificação faz aparecer 
uma nova extremidade. 

PROP. 61 — Num grafo o número total de 
malhas m excede o número de cordas de uma 
unidade 

m=c+1 (6) 


visto que existe uma malha exterior e tantas 
malhas interiores quantas as cordas. 

PROP. 62 — O número total de ramos r de um 
grafo contém o número de hastes hk e o número 
de cordas c (segundo as definições de haste e de 
cordas) 


r=h+c (7) 
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PROP. 63 — Um grafo conexo com n nós e r 
ramos tem o número de cordas c tal que 


=r-n+1. (8) 


dado que a escolha de uma árvore num grafo 
define de maneira unívoca um sistema de 
=r—-h==r—(n—1) anéis-base. 
PROP. 64 — Entre os elementos de um grafo 
G e do seu dual G' existem as seguintes relações: 


n'=m+1 (9) 
'=n-l1 (10) 
= (11) 


PROP. 65 — Entre os elementos de um grafo 
G existe a mesma relação que se verifica entre 
os elementos do seu dual G' (pois também é um 
grafo) 

=r—-n+41 
c'í=t'—-n'+1 


(12) 
(13) 


PROP. 66— A soma dos nós de dois grafos 
duais é superior em duas unidades ao número de 
ramos comum aos dois grafos 


nt+tn'=r+2=r'+2 (14) 


4. MATRIZES DE INCIDÊNCIA 


DEF. 46 — Grafo associado a um sistema físico 
linear (mecânico, eléctrico, magnético, térmico, 
etc.) é o grafo relativo às ligações entre os seus elemen- 
tos físicos. 

PROP. 67 — Um grafo associado a um sistema 
físico linear não contém informação sobre os 
respectivos elementos físicos (como consequência 
da definição de grafo). 

DEF. 47 — Matriz de incidência de um grafo 
conexo sem singularidades é uma matriz constituída 
apenas por elementos de valor + 1, —1 e O relativos 
à existência e orientação de uns elementos do grafo 
(por exemplo, ramos) em relação a outros (no exem- 
plo, nós). 

PROP. 68 — As matrizes de incidência dão 
esclarecimentos sobre a topologia dos sistemas 
respectivos mas nada dizem sobre a sua estru- 
tura física (por exemplo, fontes, impedâncias, 
etc.), pois contém a mesma informação que os 
grafos (formam-se a partir deles). 
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PROP. 69 — As matrizes de incidência corres- 
pondem aos processos de cálculo dos sistemas 
lineares evitando a formação dos corresponden- 
tes sistemas de equações. 


4.1 — Matriz de incidência dos nós 


DEF. 48 — Matriz de incidência dos nós A de 
um grafo consiste na matriz cujos elementos são defi- 
nidos por linhas em correspondência com os nós do grafo 
e por colunas com os ramos. 

DEF. 49 — Formação da matriz de incidência 
dos nós: ramos incidentes num nó determinam elemen- 
tos de valor + 1 se convergem paraonóe—1 se 
divergem do nó e ramos não incidentes determinam ele- 
mentos nulos [2]. 

PROP. 70 — A matriz À tem tantas linhas quan- 
tos os nós do grafo. 

PROP. 71 — A partir da matriz À pode recons- 
truir-se o grafo respectivo. 

PROP. 72 — Diferentes grafos com a mesma 
matriz À correspondem a sistemas lineares com 
o mesmo sistema de equações. 

PROP. 73 — Numa matriz À aparece em cada 
coluna uma vez 41 e outra vez — 1 sendo 
todos os outros elementos nulos, pois cada ramo 
(relativo a cada coluna) é convergente para um 
seu nó extremidade e divergente para o outro. 

PROP. 74 — Cada linha da matriz À é igual 
à soma negativa das restantes (devido à proprie- 
dade anterior) e portanto as linhas não são inde- 
pendentes entre si, isto é, a matriz À é singular. 

DEF. 50 — Matriz de incidência dos nós redu- 
zida A, é a matriz de incidência dos nós À sem a linha k 
(em geral a última); neste sentido diz-se que 
A é a matriz completa. 

PROP. 75 — A matriz À, corresponde a refe- 
renciar todos os elementos ão nó k (ou seja, 
equivale a tomar para as grandezas de k o valor 
de referência nulo). 


4.2 — Matriz de incidência das malhas 


DEF. 51 — Matriz de incidências das malhas O 
é aquela em que as linhas correspondem a malhas e as 
colunas a ramos do respectivo grafo. 

DEF. 52 — Formação da matriz de incidência 
das malhas: o sentido positivo das malhas é arbitrado 
no sentido directo (contrário ao movimento dos pon- 
teiros do relógio); se o sentido do ramo corresponde 
ao sentido positivo de circulação da malha o elemento 
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